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The  Finite  Element  Method  with  Iterations  for  Calculating  the 
Shapes  and  the  Frequencies  of  the  Free  Oscillations  of 
Naturally  Twisted  Propeller  Blades 

by  Z.  Ye.  Shnurov 


SUMMARY 

The  Method  of  Calculating  the  frequencies  and  the  forms  of  the 
free  oscillations  of  a twisted  rotating  propeller  blade  with  the 
various  methods  of  attaching  a blade,  encountered  in  real  designs, 
was  worked  out  by  the  Finite  Element  Method  with  Iterations.  Calcu 
lation  by  this  method  gives  quick  agreement,  has  high  resolving 
power  and  makes  it  possible  to  obtain  the  solution  in  cases  when 
it  would  not  be  possible  to  obtain  it  by  other  methods  or  it  would 
be  accomplished  with  insufficient  accuracy. 


Introduction 

The  emergence  of  vertical  takeoff  aircraft  with  rotating  propellers 
and  of  helicopters  with  rotors  having  rigid  attachment  of  the  blades  to 
hub  has  once  agiin  concentrated  attention  on  the  questions  of  the  oscilla- 
tions of  twisted  blades. 


The  increase  in  the  flying  speed  of  helicopters  and  the  necessity  of 
the  operation  of  propellers  (including  aircraft  propellers)  under  condi- 
tions very  remote  from  axial  flow  requires  the  creation  of  such  practical 
methods  of  calculating  the  frequencies  of  natural  oscillations  and  stress 


4. 


es  in  these  blades,  which  make  it  possible  with  a greater  degree  of  accur- 
acy to  determine  magnitudes  which  are  of  interest  to  the  designers.  From 
the  early  works  it  is  necessary  to  note  [1,2].  In  work  [l]  an  approximate 
method  of  calculation  was  developed,  in  which  a twisted  blade  was  simula- 
ted by  a twisted  rod:  the  plane  of  its  greatest  rigidity  was  located  at  a 
certain  angle  to  the  plane  of  rotation,  which  corresponds,  as  a rule,  to 
an  angle  on  a relative  radius  r=0.75.  In  work  [2]  based  on  the  theory 
of  thin  rods  of  Kirchhoff  and  Clebsch  accurate  differential  equations  of 
the  oscillations  of  blades  were  obtained  and  approximate  methods  of  sol- 
ving these  equations  for  a nonrotating  blade  were  demonstrated. 

The  methods  of  calculating  developed  in  works[3,4]  have  high  accur- 
acy in  solving  the  blade  model  in  question,  however,  its  twistedness  is 
not  taken  into  consideration,  as  a result  of  which  the  connection  of  the 
oscillations  in  the  planes  of  the  greatest  and  least:  rigidity  is  absent. 
These  works  also  cannot  be  used  in  the  calculations  of  plane  blades  in 
the  case  when  Inclination  of  the  plane  of  the  greatest  rigidity  in  the 
plane  of  rotation  occurs,  since  the  connection  arising  in  this  case  be- 
tween the  oscillations  in  the  planes  of  the  greatest  and  least  rigidity 
is  not  taken  into  account.  It  is  necessary  to  note,  that  in  certain  cases 
this  connection  (attachment  of  the  blade  to  the  hub  with  the  aid  of  a hor- 
izontal and  vertical  hinge)  has  an  insignificant  effect  on  the  frequen- 
cies of  the  natural  oscillations,  in  particular,  in  those  cases  when  the 
oscillations  in  the  plane  of  the  least  rigidity  predominate  in  the  actual 
form  of  oscillations. 


Of  the  works,  dedicated  to  the  calculation  of  bending  oscillations 


of  slightly  twisted  baldes,  it  is  necessary  to  not  work  [5],  in  which  & 
single-parameter  integral  equation  was  developed  for  determining  the 
forms  and  frequencies  of  bending  oscillations  of  a rotating  blade,  which 
ensures  convergence  of  the  method  of  successive  approximations  at  any 
angular  velocities  of  rotation.  But  the  case  was  examined  in  it,  when 
one  of  the  main  magnitudes  of  rigidity  becomes  infinite  and,  as  a conse<* 
quence  of  the  slight  twistedness,  the  connection  of  the  oscillations  in 
the  planes  of  the  greatest  and  least  rigidity  disappears.  An  equation  was 
obtained  in  the  work,  which  takes  into  consideration  the  twistedness  of 
the  blade,  but  the  programmed  equation  did  not  take  this  circumstance  into 
account . 

The  works  developed  in  connection  with  the  calculation  of  the  fre- 
quencies of  the  natural  oscillations  of  compressor  and  turbine  blades  can 
be  used  for  calculating  propeller  blades,  however,  the  different  assump- 
tions in  deriving  the  equations,  completely  acceptable  for  the  purposes 
of  designing  short  and  broad  turbine  blades,  do  not  make  it  possible  to 
consider  then  sufficiently  satisfactory  for  designing  propeller  blades. 
Thus,  for  example,  in  work  [6]  it  was  assumed,  that  the  curvature  in  the 
direction  of  the  greatest  rigidity  of  the  cross-section  was  equal  to  zero. 
This  assumption  is  completely  satisfactory  for  those  cases,  when  for  a 
plane  blade,  which  differs  from  the  one  being  investigated  only  by  the 
fact,  thatthe  twistedness  of  the  cross  sections  relative  to  each  other 
is  absent  and  the  frequencies  of  the  natural  oscillations  are  greatly 
different.  Such  is  the  situation  for  turbine  blades.  For  propeller 
blades  (for  example,  see  Fig.  11)  these  frequencies  are  sufficiently  dote 
and  a similar  assumption  would  lead  to  a noticeable  error. 
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In  work  [7]  a system  of  differential  equations  was  obtained,  which 
describe  the  oscillations  of  a blade  taking  into  consideration  the  con- 
nections, which  arise  with  simultaneous  bending  in  the  planes  of  the  great- 
est and  least  rigidity  and  torsion.  However,  the  method  of  initial  para- 
meters employed  in  this  work  in  the  presence  of  a field  of  centrifugal 
forces  requires  very  high  accuracy  and  does  not  give  certainty  that  the 
system  can  be  solved  at  all  possible  relationships  of  the  lowest  frequen- 
ciesof  natural  oscillations  and  angular  velocity  of  rotation  of  the 
rotor. 

It  is  necessary  to  mention  one  more  important  feature,  which  should 

be  taken  into  account  in  developing  (two  sentences  and  part  of 

another  are  illegible) 

with  aid  of  one  so-called  horizontal  hinge,  and  its  axis  cannot  coincide 
with  the  plane  of  the  greater  rigidity  of  the  blade  • bothr  at  the  site  of 
the  location  of  the  hinge  and  at  any  other  cross  eection.  This  noncoinci- 
dence leads  to  additional  connections  in  the  oscillations  of  a blade  in 
the  direction  of  its  least  and  greatest  rigidity.  Works  [1-7]  do  not 
take  such  a possibility  into  consideration. 

In  the  present  work  a method  has  been  worked  out,  which  makes  it 
possible  to  calculate,  for  a naturally  twisted  blade,  the  frequencies 
and  forms  of  the  natural  oscillations  simultaneously  in  the  planes  of  the 
greatest  and  least  rigidity  with  various  combination  methods  of  attach- 


ment to  the  hub. 
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Derivation  of  the  Equations 

The  method  expounded  below  was  obtained  by  employing  the  method  of 
three  moments  (well  known  in  strength  of  materials),  used  earlier  for 
calculating  the  deformations  of  a blade  in  one  plane  in  a field  of  centri- 
fugal  forces  [8]  and  revised  by  A.  V.  Nekrasov  [3]  for  calculating  the 
frequencies  and  forms  of  natural  oscillations  also  in  a field  of  centri- 
fugal forces  of  a plane  untwisted  blade,  performing  oscillations  only  in 
the  plane  of  the  greatest  and  least  rigidity. 


Fig.  1 


In  accordance  with  the  method  the  blade  is  represented  as  a beam, 
in  which  the  length  greatly  exceeds  the  dimensions  in  the  two  other  dir- 
ections; it  is  assumed,  that  the  points  of  intersection  of  the  main  cen- 
tral axes  of  rigidity  of  the  cross  section  lie  one  one  straight  line  and 
the  connection  between  bending  and  torsion  is  absent.  The  blade  in 
accordance  with  work  [3]  is  represented  in  the  from  of  a weightless 
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f beam,  in  each  cross  section  of  which  are  known  the  directions  of  the 

main  central  axes  of  inertia  and  accordingly  of  the  greatest  values  of 
bending  rigidity.  The  blade  is  broken  down  into  sections.  At  the  junc- 
tion of  the  neighboring  sections  are  located  concentrated  masses,  which 
correspond  in  magnitude  to  the  mass  of  the  halves  of  the  neighboring 
sections  adjacent  to  the  given  point.  The  magnitude  of  the  rigidities 
and  the  angles  of  inclination  of  the  planes  of  the  greatest  and  least 
rigidity  within  the  limits  of  one  section  are  constant  (Fig.  1).  Thus, 
the  centrifugal  force  being  applied  to  the  concentrated  masses,  remains 
constant  within  the  limits  of  one  section  and  changes  abruptly  at  the 
boundary  with  the  neighboring  section.  A zero  mass  is  selected  from  the 
conditions  determined  by  the  attachment  of  the  blade  to  the  hub.  For 
obtaining  a fixed  (jammed)  or  hinge-wise  supported  blade  the  zero  mass 
is  selected  necessarily  large.  If  the  blade  is  attached  to  an  absolutely 
rigid  real  body,  then  the  zero  mass  can  be  selected  in  accordance  with 
the  mass  of  this  real  body.  If  the  blade  is  attached  to  an  elastic  real 
body,  the  frequencies  of  the  natural  oscillations  of  which  have  the  same 
order  as  the  frequencies  of  the  natural  oscillations  of  the  blade,  then 
the  zero  mass  can  be  determined  from  the  condition  of  equality  to  zero 

of  the  sum  of  the  dynamic  rigidities  of  the  blade  in  the  corresponding 

i 

direction  at  the  zero  point  (Fig.  1)  and  of  the  dynamic  rigidity  in  this 
same  direction  of  the  attached  body.  Let  us  point  out  here,  that  for  one 
and  the  same  value  of  frequency  the  values  of  the  zero  mass  in  the 

. examination  of  the  oscillations  of  the  blade  in  plane  xOz  and  in  plane 

yOz  will  be  different. 


9. 

Let  us  examine  the  two  neighboring  sections  of  the  blade  i-j  and 
j-k  (Fig.  2).  Here  is  the  angle  of  inclination  of  the  plane  of  the 

greatest  rigidity  towards  the  plane  of  rotation  in  section  i-j;  *y»  is 
the  angle  of  inclination  of  the  plane  of  the  greatest  rigidity  towards  the 
plane  of  rotation  in  section  j-k. 

Plane  xOz  coincides  with  the  plane  of  rotation;  planes  Oz; 

02 i Oz;  Oz  are  respectively  the  planes  of  the  greatest  and 
the  least  rigidities  in  sections  i-j  and  j-k.  The  concentrated  masses 
are  located  at  points  (ij)  and  (jk).  Thus,  only  one  mass  is  located  with- 
in the  limits  of  one  section. 


Fig.  2 


Let  us  explain  more  precisely,  that  the  coordinate  system  xyz 
rotates  together  with  the  rotor,  in  this  case  axis  y is  directed  along 
the  rotor  shaft  in  the  direction  of  the  lift,  axis  z is  directed  along 
the  blade  axis,  and  axis  x is  perpendicular  to  the  first  two  and  is 


oriented  in  the  direction  of  the  rotation  of  the  rotor.  This  1b  a 


r 


10. 

so-called  general  system  of  coordinates.  Moreover,  the  so-called  local 
coordinate  system  Is  located  In  the  Investigation;  the  axes 

and  are  located  In  the  planes  of  the  cross  sections  of  the  blade 

perpendicular  to  axis  z,  are  are  directed  respectively  along  the  main 
central  axes  of  inertia  of  the  cross  sections. 

Let  us  write  out  the  equations  of  equilibrium  for  sections  i-j  and 
j-k  in  the  general  coordinate  system  xyz: 

~ — Ajy  ({-'Vy  — I (la) 

Myl.  = •»/,„-  Qx h V, ((/,;■  I 

• * i . . j w .*•  . , ^*yfl  ^ A’/A  ((  j / ( , ;)• 
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1 

1 

( 


> 
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Designated  In  the  equations  are: 

U^;  Ux  is  the  displacement  in  the  direction  of  the  axes  of  the  general 
coordinate  system; 

O,/,  Ui  is  the  oisplacement  in  the  direction  of  the  local  coordinate 
system; 

My;  Mx  are  the  flexural  moments,  the  vectors  of  which  are  respectively 
parallel  to  axes  y and  x; 

Qy;  Qx  are  shearing  forces,  the  vectors  of  which  are  also  parallel  to 
the  axes  of  the  general  coordinate  system  y and  x; 

N is  centrifugal  force,  the  vector  of  which  is  parallel  to  axis  z. 
Having  divided  each  of  the  expressions  (la)  and  (lb)  respectively 
by  N^jl^j  and  Njfcljk*  introducing  the  designations 


I 


/v~  ~fTr 

A<ri  / 


//*■ 


1 • 


(5) 


we  will  obtain 


fit  M.riJ  flj  Mm1,  — bij  (L'y,/  Uyl,)  -f-  = 0, 

fj*  M'S> -ft*  M*J,  - b,AUyj>  - U,jj)+  ^~  = 0, 


fi)  M,„  b ^L’,  , U mH\  H ^7“^  0. 

Jr-  r -•  “ ./« * - i‘.  j 

Subtracting  in  pairs  (7)  from  (6),  we  will  obtain  after  certain 


7» 


(6) 


transformations 
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t * 
S. 


•*,U .L\ 


(8a) 

(8b) 


It 


•I  • 


/.  t-: 


(9) 


Then  we  will  examine  the  deformations  of  sections  i-j  and  j-k.  We 
wiil  seek  the  deformations  in  the  planes  of  the  least  and  greatest 
rigidity.  The  differential  equations  of  equilibrium  in  the  sections 
respectively  in  the  planes  'tJOz  and  are: 

(ei:  v\y  - (A •(/;>  '= o;  <£/, u\r  - (.vur  - o.  (io> 

It  is  possible  to  write  the  solution  of  the  equations  of  (10)  in 
the  following  manner: 


Mi  /4;  sh:*.  r -fliCh;.:*:  M = /. 
where  (with  the  retaining  of  double  indexing) 


(ID 


Avj 


Ay* 


Mg i Mw  t M-ij  .. i,u 

sh  «?</  ~ Ih  au/  ’ ~ ihx,M  ~ tha.jy  ’ 

Bu/—  Mm ; B^j  = M,u  ; 

My»  _ Myj  m . _ My*  M^jt 

Shay*  ' thay*  ’ v*  shxy*  th«y»  ’ 

By*  — Mtjj ; B,j * — M,jf, 

Ml -C/iU;-.  M s-ElyUi; 


P>.yy*>lyy*>;  *vj  y»>  — pyy  y*)  /y  y*>- 


(12) 


(13) 

(14) 


Haing  substituted  expression  (12)  in  (11)  and  having  twice  integra- 
ted the  corresponding  expressions,  we  will  obtain  for  section  i-j  the 
equations  connecting  the  deformations  in  the  plane  of  the  least  rigidity 


13, 


with  the  flexural  moment,  acting  in  the  same  plane: 


or 


bu  ~ Uv, ) — dfj  Muj  + t[j  Alui  + Pj, 
bu  {U'U  — U'M ) • — t)j  Muj  — Mut  + P?j . 


(15) 


The  expression  for  section  j-k  in  the  plane  of  the  least  rigidity 
and  the  corresponding  expressions  for  the  plane  of  the  greatest  rigidity 
are  analogous.  For  example,  for  section  j-k 

(U'-J*  — Uljj)  “ 4jk  ; I 

• 1 (16> 


in  this  case 


bj* (L-'*  ~ Vi»  >~  ~/)k  Mrjt  - d)k  Mjf  + ft 


*tl’j  « -i  /-W4U 
— ; 

- t '■  • i.  - ! • 

■ ■ ' >ii  ;;  • . • 


(17) 


(18) 


is  the  tangent  of  the  angle  between  the  projection  of  the  tangent 
to  the  elastic  line  in  plane  Oz  and  axis  Oz  with  its  apex  at  point 

ij; 

: is  the  tangent  of  the  angle  between  the  projection  of  the  tangent 
to  the  elastic  line  in  plane  ■ Oz  and  axis  Oz  with  its  apex  at  point 

ii; 

is  the  tangent  of  the  angle  between  the  projection  of  the  tangent  to 
the  elastic  line  in  plane  : Oz  and  axis  Oz  with  its  apex  at  point  jk. 
Analogously  fT\:  and  ■>’,  are  tangents  of  the  angles 


14* 


between  the  projection  of  the  tangent  to  the  elastic  line  in  the  plane 
'...O:  or  \jtOz  and  axis  Oz  with  their  apices  at  points  i j ; ii;  jk; 
and  JJ  or  in  other  words,  the  components  of  the  first  derivative  of  the 
displacements  in  the  direction  of  axes  *y|  and  ( . 

Since  in  a general  coordinate  system  the  first  derivative  changes 
monotonically,  i.e.,  it  does  not  have  discontinuities,  then  for  sections 
ij  and  jk  it  is  possible  to  write 


*V 


• ’av 

*,r 


(19) 


where  $sx  and  jP  are  respectively  components  of  the  first  derivative 
of  displacements  in  the  direction  of  axes  x and  y. 


Let  us  employ  the  relation 

(20) 

It  follows  from  (19),  that 

ft,  cos  fu  ~ ft,  *»««  f„  - ft,, cot  - V)j  sin  (21a) 

ft,  sln  9„  + ft,  cot  t„  ~ ftjjtin  cos  9,».  (21b) 

We  will  obtain  from  (15)  and  (16) 

ft,  = “ 0%u ) + ft,  Myj  + Mm  ; 

ft,  - fttVs*  -Urj,)-  Mi„; 

ft,  ■“  (Uu,  — Out  ) + ft,  My  + t&tJ  Mrii  ; 
ft,,  " ~ Uw)-*Jk  Af„»  - PJk  M,j,. 

Having  substituted  in  turn  (22)  in  (21a)  and  (21b)  and  having 


carried  out  grouping,  we  will  obtain 


for  plane  xOz 


|^i)  (U*  — Um  ) cos  <fy  — b,j  — U,ji ) sin  <f,j)  •+■  [(l)j -f  d'tJ  A1,iU ) cos  — (23) 

— (0/  Muj  + d’jMm ) sin  ?yyJ— \b,k  (U^t-Uijj)  cos  <?,*  - b)k(i!,Jk  - U,jj  ) sin  * /4]-r 
-H- (djt  M.jm  4-  l)k  Mrji ) cos  •+•  (d'fk  M\jk  -j-  I’jt  Mi //)  sin  ?tk]\ 


for  plane  yOz 

I b ,j  (Uy,j  — t/:,< ) sin  -r  blf  (U^,j  — U,iit ) cos  r,/l  + |(/}y  A1,./  4”  4f> ) sin  4- 
+ (C/  *'/«/  ■+■  d'lj  Miit ) cos  ?//)  =* \blk (Ui/k  — Lh/j)  sin  -f- 
-f  & ,►  it\  - U,y/) cos  9Jt)  + \-(d)k  Mrjt  + l‘t  M jj) sin  c„  - (24) 

— \'f;t  M-jk  -f  r-k  cos  r , jt I 

Having  used  Che  relationships 

l ' rn  4 \ : o#  i - - 1\  sin  V.  {•*,  C’:sift  t — ( . ;o*  (25) 

we  will  obtain  ...(illegible)...  transformation  for  the  corresponding 
sections: 


for  plane  xOz 


(•  •• 


~\T;t  M..~- *;:i  r.  . 
.1/  • Hon*,,.  — .lis  t --  / .V;.. >#!«?,*: 


(26a) 


and  analogously  for  plane  yOz 


btfVy,  — (b,t  -j-  /'.»)  + bfkUr„ *=  (/;y  Al.i/  -f  )#in  r,>  + 

* (C/  '*•■/  4*  d,f  Mi,i ) cos  9(/  -i  (</*,  -f  l/t  sin  <i/t 
”+  (d’jt  Mi/k  -j-  l'Jk  My,)  cos  <ryi. 


(261 


16. 

Expression  (26a)  connects  the  dis- 
placements in  plane  xOz  with  the  flex- 
ural moments  in  the  planes  of  the  great- 
est and  least  rigidity  in  the  four  nei- 
boring  sections  (ii,  ij,  jj,  jk) . Ex- 
pression (26b)  connects  the  displace- 
ments in  plane  yOz  with  the  same  flex- 
ural moments  in  the  planes  of  the  least 
same  four  neighboring  sections  (ii,  ij,  jj. 

These  expressions  are  awkward  because  they  contain  values  of  flex- 
ural moments  in  four  sections,  whereas  the  displacements  are  given  in 
three  sections. 

However,  it  is  possible  to  reduce  the  right  sides  of  (26a)  and 
(26b)  to  expressions  containing  flexural  moments  only  in  three  sections. 

Let  us  note  that  moments  Mu/,  My/,  Alyy  and  M^/  on  the 

boundary  of  the  neighboring  sections  are  connected  with  each  other  by 

the  following  relationships  (Fig.  3). 

Afi/y  — Mu/  cos  (f/t  - f,y)  - M¥J  sin  (fy,  — f0); 

Mu/  ~ Mu/  Sin(?jk  — f|j)  -4-  Myy  COS  (?,y  — fy), 

whence 


!•' 

, 

f 


| (27a) 

¥ 


and  greatest  plane  in  the 


jk). 


Mu/  - Mu,  cos  (f,»  - »y)  4-  M.J/  sin  (ffk  - tyk  1 

M¥/  = — Mt//sin(i/t  — h/)  + Mu/Cot(f/t  — fls),  J 


(27b) 
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Introducing  the  designation 

A®/  = T/»  -*  »i/ 

and  having  noted  that  equally 

Mx/k  - Afu*  cot  As.,  + sin  A*,;  -=  - iW**  sin  A?,  -f  A1^t  cos  A?,.  <28) 

let  us  substitute  (27)  and  (28)  in  (26a)  and  we  will  obtain 


b‘i  U*t  — (b>i  ~ bjt)  uzj  ~ btt  d)/  cosr/y  + (l)j  cos  A?;  cos  - 

- /'y  sin  Ay,  sin  <?(y  + /)*  coss/#)  + Mr,**(rfy4  co$A9,cos«pyt-  </^sinA?Jtsin?i4)- 
- Muid?,  sin  <?/y  — sin  A?/  cos  *jy  -j-  /?y  cos  A*y  sin  f,y  -j* 

+ /Ji  sin  ?/*>—  Mat  {d},  sin  Af»  cos  ?y4  + cos  A?*  sfn  ?<»):  (29a' 

Having  accomplished  the  analogous  transformations  for  (26b),  we 
will  obtain 


M',.- -</>„ 

•;  /J-  sill  A*;CO 
V .U;ii  J't/  COf  . 


"b  7,:'|  ( , ^ ^ l/|‘y  sin*  !p^y  -l/r  ■ I/' 

V - ’ sin s ,i-  .'•! 

1 — >::•  A.'  fin  / ••■■>«  3 - 

V 

{—  </’ . sin  A.-,  sin  c;t  -f-  if;  ci'S  A^  Co- 


:oiy(.  s i 


(29b) 


Thus,  the  expressions  are  obtained,  which  connect  the  displace- 
ments of  the  three  neighboring  points  with  the  flexural  moments  in  the 
planes  of  the  least  and  greatest  rigidity  in  the  three  neighboring 
cross  sections. 


The  first  and  the  last  sections  are  singular  sections. 

Let  us  write  out  the  analogous  expressions  for  these  sections.  The 
equations  of  (15)  for  the  first  section  are: 

/>,  i (t'i.  n — l \ 1 1)  =*  Al.m  -t  fa  -i 

f K 


(30a) 
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/>,  JU  n L;ih  )^=  </*,  AJ.  iii  /, . Al, :j.  -f . . 


(30b) 


Let  us  multiply  the  first  equation  by  the  negative  value  of  sin  ; , 
and  the  second  by  cos  •;  . then  let  us  add  them,  having  used  expres- 
sions (20)  and  (25)  for  the  transformation,  and  we  will  obtain 

bu  ( ux  = cos  foi  + Af,  oi  d'i,  cos  Vo;  — Ah  oo/j",  sin  cbl  — 

- Ait  oi  dT‘t  sin  ?„  + ftx  9, 

and  having  substituted 

Af,m  = — Mi  n sin  Af,  + Af , n cos  Af,;  Aft  m — Aft  u cos  Aft  + Af,„  sin  A*,, 

we  will  obtain 

boi  — {Pm  i ” Ux  *)  =■  Af,  oo  l £,  CO*  fM  + Af,  a (< d cos  Af , COS  f (|  — 

— sin  Af , sin  f#l) — Mi  m sin  fti  — Aft  n (dj,  sin  Af,  cos  fM  -f 

■f"  co»  Af , sm  f oi)  -f*  Px  t-  (31a) 

In  an  analogous  manner,  the  equation  for  the  displacements  in 

plane  yOz  can  be  obtained,  if  we  multiply  (30a)  by  cos  f#l,  and  (30b) 

by  sin  f#t  : 

K {Uf  i — Uy  0)  = Af,  oo^o,  ?*i  -f  Af,  u (dj,  cos  Af , sin  f„  -+•  d$,  sin  Af , cos  f„)  4- 

-f-AftoofJi cos f#|  + Afj  n(— d^sinAfjSinf,,  + d^cosAflco*fo,)-r  (31b) 

The  equations  for  the  displacements  of  the  last  two  points  of  the 
blade  can  be  obtained  in  the  exact  sa?e  manner.  For  this  it  is  neces- 
sary only  in  contrast  to  the  first  section  to  employ  the  expressions 
of  (16)  and  an  analogous  expression  for  plane  *l|0z. 


As  a result  we  will  obtain 


bt-\. — Pm. r-i.  j-i)*— Afv#_i, ,-\d\_x  ,cos  f,-r.i—Af, fco| f *-i. j-j- 

+ Ah,  i-t.  ,_i  t Sin  fa— 1. 1 + Alt.  /—i.  a /J_|. , fa-i.  i + PmA  ' (32a) 


j 


19. 


.•  (U„  - Uy  - - M,,  J- 1.  ,-1  d\ f *ln  , - M,r  . /j_,.  f Sin  9,-,. , - 

- Mi.  _|.  , CO*  9,-,.  , - ,-1. , , CO*  9j_i.  , + (32b) 


It  Is  necessary  to  note  here,  that  the  systems  of  equations  for 
the  root  and  terminal (tip)  sections  differ  somewhat  from  the  remaining 
equations,  which  is  explained  by  the  peculiarities  of  the  extreme  sec- 
tions. 


Let  us  introduce  designations,  which  simplify  the  notation  of 
equations  (29),  (31)  and  (32): 


</;;  = d\,  cos  -it,:  d ’ * d'.  sin 
- -f,  v - sm  ? — •/•«  *»•#  t>: 

— /;•■  CoS  A-;  sin  f -p  /\Sin »,*; 

.*  I-  r 


- - ■ . >5-  s : ; - (/■.cos-: 

,;v  ' :,/i  /;,  Stii  A.,  Cos  •:  T ./:,A.  #111 -f  (*: 

7( >sin  .Vi,  slno„  -f  fy'cos  co»ry/  -r  >\  cos?,*: 

c/:(,  «.  o* *.in  7 ,|  • </’,  #ii!  A-.*  «.*os  ~ : 

"«•  “ _ Hn  *~r*  sin  <fy*  -+  •//*  cos  As,  cossj)t. 


As  a result  we  will  obtain  a system  of  (z+1)  pairs  written  in 
general  form,  containing  (z+1)  pairs  of  flexural  moments  acting  in  the 
plane  of  the  least  and  greatest  rigidity,  and  two  pairs  of  components 
of  the  first  derivative  of  the  elastic  line  in  fixed  coordinate  system 
xyz  at  the  zero  and  last  points  of  the  blade. 


The  general  expression  of  the  J-th  pair  of  equations  is: 
b,,Utt  — a,  Vt,  b/t  Utt  **  M,tn  d*,1,  -f  M,Jt if  + 

+ Mi,u  it  ~ Mui  — Mkjj  tj — Mm  d£'\ 

b,i  U,i  + d)  U,t  + — Myi  df,  + My),  if  -f 

+ Myt*  df  + Mui  df!)  + Mi,,  tf*  + Mut  4* . 


(33) 
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The  first  pair  of  equations  (the  root  section)  is 

+ Mr,  BO  lei  -f-  Af,  llrft*  — Ml  oo  H?  — Mind?  + fMt\  I 

— l„ i Uy  „ -f-  bal  Ut  i = Mr  ooloi  Al,i  i d*  + Mi  oo 4«  + M\  ii  d?+  fy  J (34a) 

The  last  pair  of  equations  (the  terminal  (tip)  section)  is 

b- 1 U*-  «-•  U„  = MV  - /M,,:  , /;L,; 

Mil  *-h  *->  * — Alt  i-ii  * tf! i-.t  — P„; 

ft.-i  Uy,  = Mv djLi,,  -f-  Mv*_i;  */?-!,*+  (34b) 

+ Afl;1_|;  j +•  Alj;  pyt. 

This  system,  written  in  matrix  form,  has  the  form  of  (35)  (see 
the  inset).  It  makes  it  possible  from  the  known  values  of  the  flexural 
moments  and  components  of  the  first  derivative  of  the  elastic  line  to 
find  the  magnitudes  of  the  displacements. 


For  solving  a problem  on  natural  oscillations  by  the  method 
of  successive  approximations  one  more  system  of  equations  is  needed, 
from  which  it  would  be  possible  to  obtain  the  values  of  the  flexural 
moments  and  components  of  the  first  derivative  at  the  ends  of  the 
blade. 


For  this  let  us  employ  system  (8a)  and  (8b). 


Let  us  write  equation  (8a)  somewhat  differently: 

(>/  L’,j  — bJt  UMjt  ~/4/  MrU  —fit  My,!  — 
~ /'*  Mxn  -f,t  My/o  + 


since 


21. 


r . .1; 

• 1 , \ 

— y i— 

.W:l(sinr,y  -i-  .W, 

. . cos'r,/K 

f M 

— / . 

W <in  a j.  1? 

• i*o?  “ 

-*/  i—  . 

cos  t ,A 

.'in  y-/4  -!-  AJ, 

,/*  Cl,S  T/*l 

(37) 


then,  having  equated  the  right  aides  of  equalities  (26a)  and  (36),  we 
will  obtain, after  some  transformations 


A ll'“  m,i) C0<  *•>  **  M'-'-  '!■>  + mo) cos  t.j  - ih  ?/t  - 

+ '«v*  “ ni)t ) cos  - M,i;  id\,  - Wl>)  sin  f # + „n  ?</ 

- Mxu  Oh  + »>,,)  sin  'ijt  - My*  (d)k  - m ) sm  -rjk  = 4"—  - 

v-7  ’ (38) 

Let  us  again  use  the  substitution  of  (27)  and  (28) , simultaneous- 
ly introducing  the  designations 


I."-'  !■"' . nr’.  <»,( 

•-  —,nUmm  Th  «y<'.>  lJ'  K,i  - PU  COSV,fi 

-*lv 

r--  l»i  j;  (».>  _ _ ■/*  ~ . i — li.\ 

-d,k  —m)k SiTSHv  "*/*•  p''  *=  Ptf  sin  f,,\ 


*0  li  (».)  i _ tl*» 

$</  “ lu  + m,y  “ ~ ,}yw  V'  9u  ~ 91/  cos  <piy; 

,l(y 

?/*'’  =*  0*°  + m/*  — ft) w 


(39) 


As  a result  we  will  obtain  for  plane  xOz 


* 
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2 A, 


Pit  + M, jj  (q$  -f  <3u  cos  A^  — ^ sin  A fy)  -f 
+ cos  Af4  — P^Jsin  Af„)  — Mm  Ptf  — 

— Mi/j  (qj*  + qV/  cos  A yy  -f  sin  A?y)  — Afu»  (PJJ  cos  A?4  -f  PjJ  #jn  Ayj  _ 


Analogous  transformations  will  give  for  plane  yOz 


/W,tu  Py  •+  AfT(/y  -j-  f//  COS  A 9y  — Sin  Atpy)  4-  / 

+ (/*/*  cos  Af,  •+■  Pfi  sin  A?J  + Aim  Pff  + 

+ M/y  (?yi  + if?  cos  A(py  — q* sin  AfJ  ■+■  Afu*  (P/2  COS  A?t  — P/»sin  Atp*)  = 

Q>/«  Qy</ 

A/(y  • 


(AOb) 


Thus,  equations  were  obyained  in  general  form,  which  connect 
with  each  other  the  flexural  moments  acting  at  the  three  neighboring 
in  the  planes  of  the  least  and  greatest  rigidity,  with  the  shearing 
forces  at  two  neighboring  sections  of  the  blade  (ij  and  jk^  and  direc- 
ted along  the  axes  of  the  fixed  coordinate  system. 


In  examining  the  special  root  and  tip  sections  of  the  blade,  it 
is  not  difficult  to  obtain  expressions  connecting  the  flexural  moments 
at  the  two  neighboring  points  with  the  shearing  force  acting  in  the 
root  and  tip  sections. 


25. 


Let  us  introduce  designations,  which  simplify  the  notation: 


•/;'  - ‘Li- - i?;;  cos  — q,)  (in 
P)  — fS1  r»»#  A-  . p"  cj-  i . • 

' </  C(”  j ; ---  cjisjr.  ± j 

IJ,j  cos  t Ayisdn  A-ft; 

7;  ”-'/>»  ■ </  y cos  A;.-  -r sin  A- 
/J*  M*  cos  xrt  -f-  Pj'i  Sin  As,: 

"*  7y*‘  r V./  cos  A*,  - v^j  sin  A»,; 
— J'?*  cos  at*  — A’,.'  sin  As*. 


As  a result  we  obtain  a system  of  (z+1)  pairs  of  euations  writ- 
ten in  the  general  form,  which  contain  (z+1)  pairs  of  flexural  moments 
acting  in  the  planes  of  the  greatest  and  least  rigidity,  and  two  pairs 
of  components  of  the  first  derivative  of  the  elastic  line  in  the 
general  coordinate  system  xyz  for  the  zero  and  the  last  point  of  the 
blade.  The  equations  for  the  first  and  the  last  sections,  which 
contain  these  components,  differ  somewhat  in  structure  from  the  re- 
maining equations  and  take  the  boundary  conditions  into  account.  These 
equations  are  written  out  below: 

the  general  expression  of  the  j-th  pair  of  equations  is 


M,#  tf)  + M,jj  q‘/  + P*k  — All u PJJ  — M \jj  fj'  — Afu»  PT « 

Q-r/«  Qj,tJ  , 

My  * 

Pj)  + M,jj  q‘j  ■+■  Pit  + Alur  Pfj‘  + Ali/y  q*/'  •+■ 

+ Alu*/**1  — ; 


(42a) 


the  equations  for  the  root  section  are 
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-+■  Mt, oo9oi  + M,  n PT*  — Af|  to  90?  — All  u Pi ,=  ; 

P>#  — *+■  n Pf*  + ^too?o?  + Afni  PT  = -vf-*'- ; 


K. 

the  equations  for  the  terminal  (tip)  section  are 


(42b) 


P«  + ^V-l-  *—l  P£~i.  t -r  M v_|,  , fJL|.  x — Afj,-!,  f_|  Pfl|  , — 
_ M.  . ...  Q«-l,  i . 

/»«-i.  *?,-i.-,  *=  - -tt- — ; 

'*1-1, 1 

fy*  ~f  At*-  *« * t "Wv-i.  * f<-i, « "f  ^u— i.  i-i  Pjli,  t — 


Since 


I,  i?i- 


Q.v<;  — ^ m,  ((-  — ‘«s  t'-t,);  Qyij  =*  — ^ mt  Uyt, 

> / 

then  these  are  differential  equations.  It  is  simple  to  convert  them 

into  algebraic  equations,  representing  the  solution  in  the  form 

( v,  (A)  — L.v  sin  pt\  M,j ,,,  a M,i  sin  pr. 

!\v  i'\ - l v <r:  .**.*:  .1/;  - .If .,f\r.pr. 

. • • * ' ' .* 

where  . J ; ; and  are  amplitudes  of  the 

value  of  the  functions.  Having  differentiated  twice  we  will  obtain 

• • 

(X».V  . -=  I/1'-'  2.  * • ; V * ? 

Such  a system  is  written  in  matrix  form  (43)  (see  the  inset).  It 
makes  it  possible  from  the  known  values  of  the  shearing  forces  to 
find  the  magnitudes  of  the  f lexuralmoments  in  the  planes  of  the 
greatest  and  least  rigidity  and  the  components  of  the  first  deriva- 
tive of  the  elastic  line  with  respect  to  the  ends  of  the  blade  in  a 
fixed  coordinate  system. 


i.  * 1 


Q.*1-1  * 

A'.-i, , • 


(42c) 
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The  General  Sequence  of  Calculation 

The  finite  element  method  with  iterations  developed  in  the  pre* 
sent  work  is  structured  to  a certain  extent  analogous  to  other  methods 
for  the  calculation  of  frequencies  and  forms  of  free  oscillations. 
However,  due  to  a number  of  special  features  it  is  advisable  to  state 
it  in  the  most  general  form. 

The  approximations  are  accomplished  in  the  following  order. 

The  shape  of  the  oscillations  for  the  calculation  of  the  first 
approximation  is  assumed  determined  by  the  functions  i/JS*  and  1/%?,. 

The  form  is  selected  in  such  a manner  so  that  l/f!}  — L/JS  ■■  I . . It  is 
simultaneously  assumed  that  the  frequency  of  the  natural  oscillations 
of  the  firstapppoximatlon  is  pT  ■1.  Then  the  shearing  forces 
located  on  the  right  side  are  determined  respectively  for  the  planes 
xOz  and  yOz  as 

* ■ * 

Qm  v — (1  + •*)  £ «*<  U,t\  Qm  — 1 £ «/ 

/ J I 

The  flexural  moments  in  the  planes  of  the  greatest  and  least 
rigidity  and  the  components  of  the  first  derivative  of  the  elastic 
line  in  the  root  (P„  and  ffJ  are  determined  from  the  solution 
of  system  (43). 

Then  the  displacements  and  1/mi  and  the  frequency  of  the 

natural  oscillations  of  the  first  approximation  are  determined  from 
system  (35) 


■■ 


i 
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Then  functions  UlY  and  i'j are  again  substituted  in  sys- 
tem (43)  and  the  cycle  of  the  second  approximation  is  repeated.  These 
cycles  are  further  repeated  until  the  specified  accuracy  is  obtained. 


The  method  of  successive  approximations  employed  here  leads  to 
the  obtaining  of  the  form  of  the  natural  oscillations  with  the  lowest 
frequency.  Let  us  call  it  the  frequency  of  the  first  tone,  having 
designated  the  final  form  and  frequency  in  the  following  manner 

The  determination  of  the  form  and  the  frequency  of  the  natural 
oscillations  of  the  second  tone  is  the  next  step.  The  process  of 
successive  approximations  is  repeated  as  previously,  with  only  one 


difference,  that  it  is  necessary  to  fulfill  the  condition  of  ortho- 
gonality of  the  form  of  the  second  tone  to  the  form  of  the  first. 
This  condition  is  written  in  the  following  manner: 
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where  LxV,i:.;  bW, n,  are  functions  of  the  displacements  of  the  k-th 
approximation  of  the  second  tone.  Hence,  the  functions  of  the  dis- 
placements are  written  as 

= U'/,\ id  — Ci\‘x,  (i,;  Cri  ,ii,  =»  t/y*  (in  — C:i'y,  «i», 

where  Oil’oi,  and  n are  functions  ot  the  displacements  obtained 
from  the  solution  of  system  (35); 


X mi  Umi\ id x,  si  + X-1  ^*(ii)3’/(i) 

/-(*>  iam 

c»  = nrr • 


Z mi  *i  <ii  4*  Zlmi  y?(D 

The  constant  C21  (the  constant  of  the  orthagonalization  ot  the 
second  tone  to  the  first)  varies  with  each  approximation. 


1 


~(*-n 
uy  on 


tyr*  <*>>  c,,  y/(ij 

a 

if  the  form  of  the  oscillations  of  the  first’ tone  is  such  that 
y«( ,)  ■ 1,  and 

* <S)  b xl  ,|t)  . } 

KV„„-c tom~TP*m-d#  1 


if  the  form  of  the  oscillations  of  the  first  tone  is  such  that 

(l)  ■ 1 • 

Generally  speaking,  the  process  of  successive  approximations 
must  be  structured  in  such  a manner  that  the  form  of  the  oscilla- 
tions determined  by  the  smallest  of  the  two  values  of  frequencies 
obtained  after  each  approximation  was  established  in  the  next  approx- 
imation. This  is  necessary  so  as  not  to  omit  any  form  of  the  oscilla- 
tions, since  cases  are  possible,  when 

oau-cff 

OPw-rt*. 

The  form  of  the  natural  oscillations  of  the  k-th  approximation 
of  the  second  tone  is  Mined  by  the  expressions 
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- <fe  - C!.*\v.„,)  I -f 

|1*’»  __  . / •**»  /'•*  » .\i)  v 

c rHih-*  (t  villi*  — /;,  ,jj,  , 


where  p\  means  that  normalization  occurred  with  respect  to  the 
function  of  the  displacements  in  plane  yOz 


I :.V. 


I'/  fll>  •*  (t-  wi.’l  - «.  ..  A;  li) /!.l  ,j|i. 

where  t , means  that  normalization  was  carried  out  with  respect  to 
the  function  of  the  displacements  in  plane  xOz. 


The  form  of  the  natural  oscillations  of  the  second  tone  obtained 


with  the  specified  accuracy  is  designated  in  this  manner: 

t-.ii  till  — A;  dll . 

I jri  till  — • 

Then  the  frequencies  and  the  forms  of  the  natural  oscillations 


of  the  next  tones  are  determined. 


The  alternation  of  operations  remains  the  same  as  in  the  cal- 
culation of  the  second  tone.  The  difference  is,  that  for  the  H-th 
tone  the  function  of  the  displacement  for  the  k-th  approximation  is 
written  as 


Z7J **,*)  - Cifi  x,  iM)  j (pi  (.V,  + •*) ; 

[MmN-*  1 

Urt\th  — 


1*1. 

<A>. 


the  constant  of  orthogonal izat ion  is 

Img  ^ Ini  ^ 

X mi  *•  <*»  + 5L*»<  Uyl/fiXUM) 

/-(»»  IM  W 

tai  (■! 

X mi *•  (*>  + X 

ImO  /•* 


r 


K 


31. 


the  frequency  of  the  natural  frequencies  is 


_ * 

Pm  '»»  “ mZu-\ 

— £ CjJir 

JHml 


»*  or 


M 


'ft (*. 

~B=y=r 


Uffm  — £ Cl&y 


*■1 


if  the  forms  of  the  oscillations  are  such  that  UMn'\)  • 1,  or 


_J  <*>  _ 

P*  u>  — - 


:\i=vzr 


.U-l 


„2(*> 

Py  (V)  = - - 


1 


•ll-V-l 

V 

»7Ti 


UytiS)  — Civil 


if  the  forms  of  the  oscillations  are  such  that  ■ 1. 


We  will  point  out  below  certain  of  the  details  in  carrying  out 


the  operations  in  determining  the  flexural  moments  and  values  of 
?t  and  • . 


Determining  the  Flexural  Moments  and  Components  of  the 
First  Derivative  of  the  Elastic  line 

..(First  two  lines  are  illegible) 

1)  rigid  fastening  with  flexure  in  two  mutually  perpendicular 
planes : 

2)  rigid  fastening  in  the  plane  of  rotation  and  hinged  support 
in  the  perpendicular  plane: 

3)  hinged  attachment  in  both  planes  (in  this  case  the  hinges 
can  be  located  at  different  distances  from  the  site  of  the  attach- 
ment of  the  blade) . 


By  analogy  with  a helicopter  rotor  we  will  call  these  hinges 
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horizontal  hinges  and  vertical  hinges. 

Naturally,  the  different  conditions  of  attachment  introduce 
differences  into  system  (43),  which  is  written  in  general  form  so 
that  it  would  be  possible  to  examine  the  case  of  simultaneous  attach** 
ment  in  two  planes  or  simultaneous  hinged  attachment  in  these  same 
planes. 

Let  us  first  examine  the  simpler  and  more  general  case  of  a 
blade  completely  fastened  at  the  root. 

The  boundary  conditions  in  the  root  ere:  ^ » 0;  ’?»  — 0 

I ml  I ml 

£/n,U,  = 0. 

/•  o * 

Generally,  system  (43)  represents  on  the  left  side  four 
three-band  submatrices,  joined  by  common  unknowns.  For  the  solution 
we  employed  the  Gauss  method  used  by  A.  V.  Nekrasov  in  [3]  for  sol- 
ving a system  representing  one  of  the  indicated  submatrices. 

The  algorithm  is  shown  below,  by  which  the  matrix  is  resolved,  with 
the  necessary  explanations  prefaced  to  it. 

Let  us  examine  system  of  equations  lx  and  ly  [see  (43)  in  the 
inset].  Since  the  blade  is  fastened  "0),  it  is  possible 

to  write 


(44) 
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where 

,mt  . | \ 
f8- -(#»+•?)£  <•! tv. 

imt 

«--p* 

(■1 

Let  us  represent  the  flexural  moments  In  the  following  manner: 


M,t  #o  **  7^  o + ^ o M,'  ii  + o Mi  u I 
Afi*>“»7ia  -ft.#  » 


(45) 


where 


| Pjl  got  + P'i'  got  . i P!l£m — Pi"  En  .. 

go i got  + £oi  got  got  got  ~r  fcigil 


f pr 

_rr,  ot\  _»» 

fTn!  — g'". 

rr-  eti  - P-“  r-  1 

> * 1 .>  • . ■■  9 . 1 

(46) 

1 

V:  ***■  r 

V 4,'vi  i.  -l 

‘:r-  r!  rr  -i  ..  ■ 1 

i.*’*  ~ » ' i 

1 

r 

,*» 

- X • 

' Fv  -L  £•;  Fi 
_ v ",y  • • 

r:  i»y  *i  L •. 

j»<;l 1\\  —r  gy.  < j 

Then  let 

us  examine  the 

system  of  equations: 

.U.  — 

: Cl*  - ..  1 ,s-  • 

Mi,.- p{';  ~ 

* : * * -4  i 

.11.  u /Vm1-  .11.  r 

J-  .1/  •• 

.1/  (V  Pft;' - .11.  !•  L>  ,11;.; 

I<47> 

where 

/•',  = - (//•-;-* 

,S»I— 

\ i k .• 

■v');  vxr  xv 

; >*,-*> -y  • : 

*'•1 

H i»iJri  + S,Jr,). 

If  we  then  substitute  (45)  in  (47)  and  carry  out  grouping, 
then  it  Is  possible  to  obtain  an  expression  analogous  to  expression 
(44); 

.‘I.,  „ /i**  -j*  Mr  s; Pa  — Mi  u gl'-  — Al» u P»‘  — ff  | 
ill.  ii  gi*"  t MKtiPfl  r Mfugi’’  + M<n  ff, . „ I 

*i  ►*'  J ; » <*> 


v * (48) 

* 1 
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where 


ti1'  - gV + Hi  Ko-  p»r\o\ 
gi*  — gV‘  4-  Pm  + /£?£.,*} 


it  ‘ ^gV1—  P'oihn  +,  ^oC*4  o'. 

* ;.-r. 

#P*  — gP  + Pi\U»  + Pit 

Fr-n-p&  T,*~>»7w. 


(49) 


These  operations  are  sequentially  repeated  for  each  pair  of  equations. 
For  the  j-th  pair  they  are  written  in  the  following  manner: 


where 


if 


if 

py- 

*V  = 


•>«/  = 


= 8?  + Pi)  Ki  — Ptj  i’,1 ; 
= g'/‘  + Pi) *v  + Pi/  i>,i ; 
■ P*J—  Pi)  Tv  + Pi)  Mi ; 

tf'gi'r  + gftf' ; 
P?  gi''-P7-j T’ . 


i 7’  "■£/' — Pi/iy—  PlJ’^u  ; 
g7'  — g7  4-  Pi)  tv  4’  Pi?  fu< 
n'^F’j-Pi'/Ms-Pihu, 

j _ pur-rfgf 

grgr+tf'gr'1 

pr'gr+p?gr 


•*.i 


g?  gt  ■ 4 -gV  gP 

+ <,r  pf 

,/s*  -»»,•  . »;•  !*•  • 

/l . /J  * — — /» , Jl.  • 


0|,  — — 


V/  f 


f/ 


•ii 


-ft*  _c»*  I _i;*  -rr*  * 
9i  ?l  + ?/  <Jl 

-rf*  cl*  -»•*  c»* 

9‘  ‘ I — ' l 

-fS*  -r»*  , jf  i.*  • 

f I *f  Qi  fl 


(50) 


(51) 


The  expressions  for  the  flexural  moments  are: 


■■  » v ' °v  -^\jj  - ,ju  Aim  ™ 71/  + £../  /H> _//  -f  Mi//‘. 

-""!3Jr,  - v~);  */— >Trrr-  . 


(52) 


As  a result  we  will  obtain  equations  connecting  in  pairs  the  two  flexural 
moments  in  the  plane  of  the  greatest  rigidity  with  the  two  flexural  moments 
in  the  plane  of  the  least  rigidity  right  up  to  ...(illegible) 


By  carrying  out  a similar  operation  for  the  z-th  pair  of  equations  . 
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(illegible) .. .we  will  obtain 


' • , , r r*  I 

M:-\.  t-1  £:-\~  *'lu—  4.  : I 


Ait  i - 1 : S?--.  -■ 


(53) 


since  at  the  free  end  of  the  blade  , ■=  Ah,-i.  • 0. 


Employing  the  discussed  algorithm,  it  is  possible  to  write,  that 

Ah,- 1,  ,.i  — ,_i  — f.,i.  i, 

where  7:*-i  and  y,._i  are  also  determined  according  to  the  formulas  of  (51). 


Having  obtained  from  (53)  the  values  of  moments  Ahz-\  and  .--i, 

it  is  then  possible  in  in  reverse  order,  solving  in  turns  the  systems  of 
two  equations  of  the  type  of  (52),  to  obtain  the  values  of  the  flexural 
moments  in  all  cross  sections  of  the  blade: 


Aitji  “ 


( rr  Ml 


ii 


py— m 


tfts'+irer  . ■ 


There  remains  the  task  of  determinign  the  values  of  ^ and 

It  is  possible  to  obtain  tern  from  the  solution  of  equations  (z+1)  (z+1)  . 

x y 

However,  these  values  are  not  needed  in  the  successive  approximations  for 
a standard  blade.  The  exception  is  the  case  when  the  end  of  the  blade  is 
not  free.  In  this  case  the  penultimate  pair  of  equations  should  be  rear* 
ranged  depending  on  these  conditions,  and  accordingly  the  moments  on  the 
end  of  the  blade  and  the  components  of  the  first  derivative,  of  the  elastic 
line,  determined  from  system  (43)  will  participate  in  the  successive  approx- 
imations. We  will  not  dwell  on  this,  since  the  corresponding  operations 
are  analogous  to  those  untertaken  with  different  boundary  conditions  in  the 
root  part  of  the  blade;  these  operations  are  discussed  below. 
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Another  case  is  rigid  attachment  in  the  plane  of  rotation  and  hinged 
attachment  in  the  vertical  plane. 


In  comparison  with  the  first  case  the  boundary  conditions  vary  some- 


what, notably 


£-0; 

mi  “ Of  3 mi  Uy,  =*  0. 

Ia0  ill 


In  consequence  of  hinged  attachment 

yW:»cos?g,  — AK.ooSin?,,  = U. 


(55) 


Hence,  the  following  order  of  the  resolution  of  system  (43)  follows: 
in  system  (43)  equation  1^  is  replaced  by  equality  (55),  and  all  the  sub 
sequence  steps  are  retained  unchanged  in  comparison  with  the  case  of  total 
attachment.  The  expressions  of  (46)  are  also  retained,  it  is  necessary  as 
soon  as  possible  to  place  in  them 

A'®*'  — cosc01:  in  »»,;  P'r  - Pf*  = FI  - 0. 

The  third  case  is  the  hinged  attachment  in  both  planes;  the  hinges 
can  be  located  at  different  distances  from  the  zero  point. 


System  (43)  was  obtained  by  assuming  the  equality  of  the  derivative 
of  the  elastic  line  on  the  right  end  of  the  i-th  section  and  on  the  left 
end  of  the  (j-k)-th  section.  The  hinge  at  the  j-th  point  violates  this 
condition  and  one  of  the  equations  of  the  (j+l)-th  system  (43)  should  be 
replaced  by  an  equation  which  takes  into  account  the  presence  of  the  hinge. 
If  the  axis  of  the  hinge  is  located  in  the  vertical  plane  the  following 
condition  should  be  fulfilled 

• ••  .m . , >*»?«-».«.  (56) 

If  the  axis  of  the  hinge  is  located  in  the  horizontal  plane  the  fol- 
lowing condition  should  be  fulfilled 


! 


! 
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.!/•  j-ir?  , - M „ o. 


Thus,  In  the  sequential  transition  from  the  first  pair  of  equations 
to  the  last  the  (j+l)-th  system  should  be  solved  (for  example,  in  the 
case  of  a vertical  hinge) : 


M’Jj  KS  + M’M  ~r  Mi/j  Ej'  + Mu,k  P?‘  Fj  : 

.M.j/  cos  zJt  — M-jj  sin  cik  — 0. 


With  this  order  of  the  resolution  of  system  (43)  all  the  expressions 
of  (51)  are  retained,  it  it  is  assumed  that 

= cos  'ijk : gl;:  = Kn  fy;  /•/'  = /•£■*  A**' « (j. 

In  the  case  of  the  positioning  of  a horizontal  hinge  at  any  i-th 


point  the  (i+l)th  system  is  solved: 

g“  —*  Mnf  gp  — Mfjj Pj*1 = Ft 

M^u  sin  f „■+  Mui  COS  1tJ  - ti, 

and  there  will  accordingly  be  obtained  in  the  expressions  of  (51) 


gT  — cos  ft,;  g?  - sin  f„;  P?  K p?  . Ff*  • Q, 

If  hinges  are  combined  at  any  point  k (including  at  the  zero  point), 


then 


*»  0;  Muk  > 0. 


For  conserving  the  algorithm  it  is  necessary  in  the  expressions  of 
(51)  to  assume 

^Sfi  ■*  " Pit-i  ■ FCn  ■ /'*♦!  ■ 0. 


The  appearance  of  hinges  imposes  the  requirement  for  one  more  opera- 
tion in  solving  system  (43)  - this  is  the  determination  of  the  component 
of  the  first  derivative  of  the  elastic  line  in  the  plane  perpendicular  to 
the  axis  of  the  hinge.  The  requirement  can  be  fulfilled  by  solving  the 
special  equation  for  the  section,  the  left  boundary  of  which  is  combined 
with  the  hinge. 


✓ 

K 
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In  accordance  with  (15)  and  (16)  we  have 

bu  (Ujjj  — (/,«)**=  di'/Mu/  -f  lij Mut  + ??<; 
bii  (^4/  — Um)  *“  dy  Af.ji  -f  ?i». 

Having  multiplied  the  first  equation  by  the  negative  value  of  $|nc., 
and  the  second  by  cos  <?,/  , let  us  add  them,  employing  (28) . We  will  obtain 


bu  Wx)  ~ uxi)  = M\U  l)j  CO*  ?jy  + Mrjj  ( d)j  cos  A?y  cos  *,y  — s»n  A=y  sin  s;y  d]y)  - 
— Mm  sin  9ij  — AUjMj  *m  *ry  cos  ?iy  + (Tij  cos  A* y sin  ?,y)  4-  'iti.  ( 57) 

From  (7)  we  will  obtain 


or 


" >1  ...  L ‘ ,)  - / , .1 1 y - /. , y,  ~ 

•'  (/ 

•*".  ' .•/:  ys!ns;  -r  .11,1/ cos?,. I — 

'a'  / 

•••••'  \ 


After  the  transformations  taking  into  account  the  substitutions  of 


(41)  let  us  finally  determine 


(58a) 


Carrying  out  analogous  operations  it  is  possible  to  find  the  expres- 


sions for  determining  pf  ; 


•*yi 


.11,, i & - MwP?  — .11,,,  - .1!i»  PT*  - 


"AT 


(58b) 


It  is  easy  to  see  that  equations  (58a)  and  (58b)  do  not  differ  in 
structure  from  analogous  equations  for  the  root  section  of  the  blade  (42b) 
For  determining  components  jiA,,  and  :y(,  it  is  sufficient  to  use  the  ex- 
pressions 

- Af.Mfat -r + /M„ PT-  + ; 

P?1— ^issfff — *-• 

If  a combined  hinge  occurs  somewhere,  then 
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P,I  = — P*f  — M\h  P?'  + . 


(59) 


Determining  the  Displacements 


System  (40)  is  used  for  determining  the  displacements;  this  system  in 
essence  represents  two  unconnected  systems,  one  for  the  plane  xOz,  the 
other  for  plane  yOz.  Taking  the  boundary  conditions  into  consideration 

im* 

For  this  the  displacements  at  each  point  are  represented  in  the  form 
of  a binomial 

UMl  + U„  -U,  , + 1/„;  (60) 

where 

In  accordance  with  this  the  determination  of  the  displacements  breaks 
down  into  two  steps. 


In  the  general  case,  in  the  first  step  the  displacements  Utl  and 

* 

Uy,  are  determined  from  simple  reccurrence  formulas: 
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( AITi/  d'j  — 


— M..i ; -m 


Mtj:  (/'+  M-kt  df  — Mm  d’j  — 

* * 

«'  • .i  fl  — OjUyl): 


(61a) 


’ f *• 


(61b) 


Then  U and  U are  determined 


l 


v r 


t •/ 

::  w. 


With  the  presence  of  hinges  in  accordance  with  (57)  the  formulas 
vary  only  for  the  {joint  following  after  a hinge. 


If  a vertical  hinge  is  located  at  point  j , then 

4*  fa  + rf*  - M fa  - M:tt  <C*  -r  l>jt  ,).  (63a) 

If  a horizontal  hinge  is  located  at  point  j,  then 


l'k  /-‘  *'**-I‘ a"  h.t  ’l  r*  (63b) 

The  subsequent  process  occurs  according  to  the  same  formulas  (61a) 


and  (61b). 


In  the  case  of  the  placing  of  a hinge  at  the  zero  point  formulas  (63a) 
and  (63b)  change  accordingly.  In  the  case  of  the  placing  of  both  hinges 
at  the  zero  polntthe  formulas  are  simplified 


Ut , - J-  (M,i.  - Af,„  + K .'5 

"01 

Uy , “ -g — (Af,,u  dil  + Miu  + Py  •)■ 
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The  possibilities  of  the  method  are  not  examined  In  the  present  sec- 
tion In  detail,  however,  the  figures  presented  below  make  it  possible  to 
give  It  a certain  evaluation. 


41. 


I 

I t 

t 

b 

1 


The  most  important  characteristics  of  the  method  are  the  required 
value  of  accuracy  for  obtaining  a satisfactory  solution  and  the  calcu- 
lation time.  It  appears  that  the  method  discussed  above  requires  calcula- 
tions with  relatively  low  accuracy,  and  the  number  of  approximations 
needed  is  such,  that  for  the  purposes  of  calculation  practically  any 
electronic  computer  can  be  used.  The  results  of  calculations  presented 
below  were  obtained  on  a M-20  computer. 

All  the  calculations  were  carried  out  with  respect  to  a hypothetical 
blade,  the  prototype  of  which  was  the  blade  of  the  tail  rotor  of  the 
Mi-6  helicopter.  The  diameter  of  the  rotor  was  taken  equal  to  6.8  m and 
its  rate  of  rotation  n = 680  r/min.  The  other  characteristics  of  the  blade 
varied  depending  on  the  purposes  of  the  calculation.  These  characteris- 
tics are  shown  in  Fig.  4.  Also  show  here  are  how  these  data  were  repre- 
sented for  calculation  on  the  computer. 

Fig.  5 shows  the  number  of  approximations  n and  the  time  t , re- 
quired for  accomplishing  these  approximations,  depending  on  the  accuracy 
, which  is  defined  as  the  greatest  difference  between  the  values 
of  the  deformations  at  any  point  of  the  blade  in  two  successive  approxi- 
mations. It  is  possible  to  see,  that  the  setting  angle  of  the  blade  with 
respect  to  the  plane  of  rotation  somewhat  affects  the  number  of  approxi- 
mations. 

Fig.  6 shows  the  change  in  the  frequency  of  the  natural  oscillations 


of  the  blade  depending  on  the  accuracy  of  the  solution. 
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It  turns  out  that  beginning  with  an  accuracy  of 
the  value  of  the  frequency  ...(rest  of  sentence  and  next  paragraph  is 
illegible) .... 

This  characteristic  becomes  especially  important  when  the  actual 
frequencies  of  the  natural  oscillations  are  close  to  each  other.  In  this 
case  the  values  of  the  natural  frequencies  obtained  by  calculation  appear 
to  be  sensitive  to  the  accuracy  of  the  definltioncf  functions  employed  in 
the  calculation. 

This  circumstance  compels  the  calculators  to  somehow  modify  the 
calculational  methods,  to  increase  the  accuracy  of  the  definition  of 
functions  and  the  accuracy  in  the  intermediate  operations.  And  this  in 
its  turn  is  frequently  limited  by  the  possibilities  of  computers  or  in 
any  case  significantly  increases  the  calculation  time. 

A certain  amount  of  uncertainty  in  the  reliability  of  the  obtained 
results  does  not  favor  the  successful  application  of  these  methods. 

yponaiMUHcn  of  t-htft 

While  not  carrying  out  a systematic  investigation  of  the| method 
examined  above,  let  us  illustrate  them  by  a comparison  of  it  with  the 
Bubnov-Galerkin  Method  in  the  form  in  which  it  was  employed  in  solving 
analogous  problems,  for  example,  in  work  [l]. 

Fig.  7 shows  the  dependence  of  the  frequency  of  the  natural  oscilla- 
tions of  a blade  fastened  to  a hub  with  the  aid  of  one  horizontal  hinge 
on  the  angle  of  Inclination  of  the  plane  of  the  greatest  rigidity  to  the 
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does  not  make  it  possible  to  find  frequencies  over  a broad  range  of  set- 
ting angle  values  (14°— 44°).  Moreover,  it  is  necessary  to  emphasize  that 

in  a whole  series  of  cases. .. (remainder  of  sentence  is  illegible) 

It  is  also  necessary  to  consider  that  the  Bubnov-Galerkin  Method  is 
favored  by  its  simplicity  and  by  the  fact  that  it  makes  it  possible  to 
employ  ready  developed  operation  elements,  from  which  a future  methodol- 
ogy can  be  compiled.  For  example,  let  us  cite  the  results  of  the  cal- 
culation of  the  dependence  of  the  frequencies  of  a plane  untwisted  blade 
attached  to  a hub  by  a horizontal  hinge  on  its  setting  angle. 

The  dependences  obtained  employing  the  Bubnov-Galerkin  Method  and 
the  method  of  successive  approximations  (Fig.  8)  very  satisfactorily 
agree  for  the  first  two  tones  of  the  natural  oscillations  of  the  tail 
rotor  of  a helicopter.  For  higher  tones  their  divergence  attains  20%. 

The  calculations  cited  for  blades  with  a different  combination  of 
parameters  make  it  possible  to  more  acs;urately  define  certain  assump- 
tions and  simplifications  established  in  blade  design  practice. 

Thus,  it  is  usually  assumed,  that  the  frequencies  and  the  forms  of 
the  natural  oscillations  of  blades  with  a typical  construction  (for 
a helicopter  rotor)  of  the  attachment  of  the  blades  to  the  hub  (horizon- 
tal and  vertical  hinge)  very  slightly  depend  on  the  design  twistedness  of 

y 

the  blade  and  the  setting  angle  relative  to  the  plane  of  rotation.  It  is  j>' 

r 

evident  from  Fig.  9 that  both  for  a twisted  as  well  as  for  a plane  blade 

i . 

the  frequency  of  the  natural  oscillations  with  preferential  deformations 
in  the  plane  of  rotation  rather  noticeably  depends  on  the  setting  angle 
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Fig.  9- 


Fig.  10, 


Fig.  9.  Key:  1)  r/min;  2)  [osc/min];  3)  tone;  4)  blade  has  design 
twist;  5)  blade  is  not  twisted. 

Fig.  10.  Key:  1)  [osc/min];  2)  tone. 
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Fig.  11.  Key:  1)  [osc/min];  2) 
taking  the  design  twist  into 
account;  3)  without  taking  design 
twist  into  account;  4)  9-th 
harmonic;  5)  tone;  6)  [r/min]. 
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within  the  usual  limits  for  helicopter  blades.  It  is  obvious  that  this 
sensitivity  to  setting  angle  depends  on  blade  characteristics,  in  parti- 
cular, on  the  relationship  of  the  rigidities  in  the  plane  of  the  least 
and  greatest  rigidity.  The  results,  shown  in  Fig.  10,  of  the  design  of 
a blade  with  the  greatest  rigidity  reduced  by  10  times  as  compared  with 
the  blade,  represented  in  Fig.  9,  illustrate  this  fact. 

In  conclusion,  let  us  illustrate  the  advisability  of  taking  blade 
twist  into  consideration  for  calculating  the  frequencies  and  the  forms 
of  the  natural  oscillations  with  respect  to  the  blades  of  one  of  the 
aircraft  propellers. 

Fig.  11  shows  the  resonance  diagram  for  such  blades,  obtained  by 
taking  into  account  twist  in  the  cetrifugal  force  field,  and  also  pre- 
sented here  is  a resonance  diagram  which  does  not  consider  the  design 
twist  of  the  blade,  l.e.,  because  this  is  done  in  design  practice.  One 
a.asily  be  convinced  of  the  advisability  of  calcualting  by  a methodo- 
lof  , ;hich  takes  blade  twist  into  consideration. 
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Flutter  Analysis  of  a Helicopter  Rotor  in  Flight 

by  V.  V.  Nazarov 

Summary 

Methods  of  flutter  analysis  of  a helicopter  rotor  in  flight 
have  been  developed;  these  methods  make  it  possible  to  consider 
blade  elasticity,  the  effect  of  air  compressibility  and  the  plia- 
bility of  the  cyclic  pitch  control  system.  A description  of  the 
methods,  a block  diagram  of  the  program  and  certain  results  of 
flutter  analysis  of  rotors  are  presented. 


.V  - the  number  of  rotor 
blades; 

«*  - rotational  speed  of 
the  rotor  [l/s]; 

{ (»,)  - moment  of  control  rod 
relative  to  the  fea- 
thering (flapping) 
r hinge; 

* m -if-  m gear  ratio  of  the 

blade-flapping  con- 
trol; 

t„  k.  - coefficients  of  the 

total  rigidity  of  the 
of  the  elements  of  the 
check  respectively  of 
the  pitch  and  lateral 
control  and  the  col- 
lective pitch  control 
from  the  control  stick 
to  the  cyclic  pitch 
control 


- dimensionless  flying  speed; 

- rigidity  of  blade  in  flex- 
ure; 

b - blade  chord; 

GT  - rigidity  of  blade  in  torsion; 
? *•  ty  - angles  of  Inclination  of  the 
cyclic  pitch  control  respec- 
tively during  deformation  of 
the  cables  of  the  pitch  and 
lateral  control; 

* , -complete  angle  of  turn  of  the 
transverse  cross  section  of 
the  j-th  blade  due  to  the 
kinematics  of  the  blade-flap- 
ping control  and  the  elastic 
deformations  of  the  blade; 

- ...(Illegible)... 

- ...(Illegible)... 

.v,  — --  - relative  position  of  the 

axis  of  the  feathering  hinge. 


Conventional  Designations 


M „« 
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- coefficient  of  total 
rigidity  of  elements 
of  cables  from  the 
blade  to  cyclic  pitch 
control: 

A 

• - relative  distance  from 

* transverse  cross  sec- 
tion of  the  blade  to 
the  rotor  axis; 

- azimuth  of  the  control 
rod  of  the  j-th  blade 
of  the  rotor; 

W - flying  speed  of  the 
helicopter  [m/s]; 

j - rotor  blade  number 

.'•I,  — - circular  Mach  number 

of  rotor; 

/j  - azimuth  of  the  j-th 
blade; 


The  broken  line  designates  the 
derivative  with  respect  to  the  radius, 
the  dotted  line  designates  the  deriva- 
tive with  respect  to  time. 


Introduction 

The  purpose  of  this  work  is  the  creation  of  a method  of  analyzing 
a rotor  for  flutter  in  flight,  which  makes  it  possible  to  consider  both 
the  elasticity  of  the  rotor  blades,  as  well  as  the  pliability  of  the 
control  system,  kinematically  connecting  the  motion  of  all  the  rotor 
blades . 


This  refinement  of  the  analytical  scheme,  in  comparison  with  the 
analytical  methods  [3]  developed  earlier  by  other  authors  and  in  com- 
parison with  the  analytical  method  discussed  in  [7]  was  evoked  by  the 
following  circumstances.  First,  the  experiments  with  dynamically  simi- 
lar models  of  rotors  of  large  diameter  showed,  that  the  effective  cen- 
tering of  the  blades  (see  [3])  for  such  rotors  is  not  a comprehensive 
characteristic  of  distributed  antiflutter  balancing.  The  balancer,  which 
changes  the  effective  centering  ly  one  and  the  same  value,  is  mounted 
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at  different  cross  sections  of  the  blade;  it  cahnges  the  critical  num- 
ber of  revolutions  of  rotor  flutter  differently.  This  can  be  explained 
only  by  the  effect  of  blade  elasticity,  which,  consequently,  should  be 
taken  into  account  in  flutter  analysis.  Secondly,  the  pliability  of 
the  control  system  cables,  which  is  usually  taken  into  account  by 
introducing  an  isolated  blade  of  effective  rigidity  into  the  calculation, 
cannot  be  taken  into  account  in  the  case,  when  the  lateral  control  of 
the  helicopter  is  accomplished  not  by  means  of  the  cyclic  pitch  control, 
which  is  characteristic  for  helicopters  of  transverse  arrangement.  In 
this  case  the  rigidity  of  the  attachment  of  the  cyclic  pitch  control  with 
respect  to  the  lateral  inclinations  is  extremely  great,  and  even  in 
the  hovering  mode  the  analysis  of  rotor  flutter  does  not  agree  with 
the  calculation  of  the  isolated  blade  [6].  Moreover,  as  is  shown  in 
[7],  the  calculational  set  up  , which  makes  it  possible  to  calculate 
rotors  of  high-speed  helicopters  and  rotary-wing  aircraft  of  different 
arrangements  in  different  modes,  for  example,  in  the  mode  of  large  , 
it  is  necessary  to  include  in  it  the  calculation  of  the  effect  of  the 
compressibility  of  air  on  the  aerodynamic  characteristics  of  the  air- 
foils and  to  make  it  possible  to  take  into  account  their  periodic 
change  under  flight  conditions  more  accurately,  than  by  simple  aver- 
aging. This  becomes  necessary  in  calculations  with  large  numbers. 
Moreover,  the  purpose  of  this  work  was  the  achievement  of  a method  with 
all  its  variations  in  a single  algorithm  and  in  a single  program  for 
a computer,  of  sufficiently  fast-response  for  carrying  out  parametric 
investigations . 


Assumptions.  General  Equations 
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In  the  investigation  of  flutter  it  is  necessary  to  take  the  fol- 
lowing factors  into  consideration  as  a totality: 

- the  elasticity  of  the  blades  and  the  rotor  control  system, 

- the  high  subsonic  tip  velocities  of  the  blades, 

- the  large  nonlinear  forces,  close  in  character  to  the  forces  of 
...(Word  illegible)...  dry  friction  arising  during  oscillations 

of  the  blades  in  the  feathering  hinge,  loaded  with  centrifugal  forces, 

- the  nonsteady-state  nature  of  the  flow  around  the  blades,  when 
disturbances  remain  for  a long  time  close  to  the  spinning  disk. 

The  following  assumptions  are  introduced  when  developing  a method 
for  analyzing  rotor  flutter. 

1)  the  hypothesis  of  plane  transverse  cross  sections  is  assumed 
valid,  in  accordance  with  which  the  blade  is  replaced  by  a beam  with 
variable  elasto-mas  linear  characteristics,  to  which  is  fastened 
nondeformed  transverse  plate-sections; 

2)  the  blade  is  considered  untwisted;  we  disregard  its  static 
flexural  and  torsional  deformations; 

3)  the  mass  of  the  blade  is  assumed  distributed  in  the  plane  of  the 
chord ; 


r 

i 
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4)  before  the  appearance  of  the  oscillations  the  plane  of  the  blade 
chord  Is  assumed  perpendicular  to  the  rotor  axis; 

5)  the  flapping  and  drag  hinges  are  assumed  to  be  located  at  one 
stationary  point,  at  which  the  boundary  conditions  at  the  blade  root 
are  assigned; 

6)  the  hypothesis  of  stationariness  is  assumed  valid  in  determining 
the  aerodynamic  effects  on  an  oscillating  blade; 

7)  the  assumption, approximately  valid  for  the  majority  of  blade 
designs  is  introduced;  this  assumption  states  that  before  the  appear- 
ance of  oscillations  the  axis  of  rigidity  of  the  blade  is  a straight 
line,  coinciding  with  the  axis  of  the  feathering  hinge. 


In  the  assumptions  enumerated,  the  differential  equations  of  the 
movement  of  a blade  element  of  unit  length  can  be  written  in  the  fol- 
lowing form  (see  [2],  [3]) 


<£/«*)'  +.'|(W)'  J mti,J  "T /I  — 0 
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where  Pr  and  are  a linear  aerodynamic  force  and  moment,  acting  on 
an  airfoil,  oscillating  in  a plane-parallel  flow. 


The  expressions  for  Pr  and  , derived  on  the  basis  of  the 
hypothesis  of  stationariness  [l]  have  the  form 
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where 

V-wr+lTccrsx;  K,  - - IP  sin  y;  i.s-g-V"  ; X'f=-  4^*-. 

. • ' **  < 

The  differential  equations  of  (1)  are  used  for  investigating  the 
oscillations  of  blades  during  flutter,  which  arise  and  develop  near 
steady-state  flywheel  motion. 


The  aerodynamic  coefficients  of  the  profiles  of  the  transverse 
cross  sections  of  a blade  v.  and  which  enter  into  the 

equations  of  (1),  essentially  depend  on  the  magnitude  of  the  local 
Mach  number,  i.e.,  on  the  normal  velocity  component  of  flow  V,  incident 
on  the  blade.  Since  the  velocity  of  the  incident  flow  V varies  periodi- 
cally with  time,  then  the  aerodynamic  coefficients  are  periodic  func- 
tions of  time  and  at  each  given  moment  they  are  different  for  different 
cross  sections  of  the  blade.  Thus,  the  equations  of  (1)  represent  a 
system  of  differential  equations  in  partial  derivatives  with  variable 
coefficients,  which  are  periodic  functions  of  time.  In  forming  the 
boundary  conditions  for  equation  system  (1)  a number  of  assumptions  is 
introduced,  which  schematize  the  rotor  control  system.  In  this  case, 
the  usual  type  of  rotor  is  examined  with  individual  hinge  attachment 
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of  each  blade  to  the  hub  and  with  cyclic  pitch  control.  It  Is  assumed, 
that  all  the  elements  of  the  control  cables  are  weightless  and  the 
elasticity  of  the  control  cable  elements  Is  taken  Into  consideration  by 
Introducing  appropriate  elastic  elements.  The  friction  forces  (see  [3]) 
are  taken  into  account  only  in  the  feathering  hinge,  since  the  magnitude 
of  the  moment  of  these  forces  Is  of  approximately  of  the  same  order  as 
the  magnitude  of  the  moment  of  the  external  forces  acting  on  the  blade. 
It  is  possible  to  show  [5],  that  with  these  assumptions  the  boundary 
conditions  in  the  root  of  each  blade  take  the  form  (see  [5]) 


where  the  magnitude  of  expresses  the  connection  between  the  rotor 
blades  through  the  cyclic  pitch  control  and  is  equal  to  the  displacement 
along  the  vertical  of  the  lower  point  of  attachment  of  the  control  rod 
of  the  J-th  blade  of  the  cyclic  pitch  control. 


It  is  possible  to  show  [5],  after  having  compiled  the  equations  of 
equilibrium  of  the  cyclic  pitch  control  under  the  effect  of  all  the 
elastic  forces,  that 
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It  is  shown  in  work  [7],  how  the  equation  of  flutter  by  the  Bubnov- 
Galerkin  Method  reduces  to  a system  of  usual  equations  with  periodic  coef- 
ficients. Taking  the  pliability  of  the  cyclic  pitch  control  into  account 
the  equation  of  flutter  for  the  j-th  blade  takes  the  form 


Cx,  + wD[t)jij  + -s  [B  (0  + C*  1 x,  + . 

+ A | x , - £ (C , 4-  20,  sin  sin  + 20,  co*  % cos  *„)  x . J = 0 

</-  0.1.2 ; A*  — 1);  (5) 

here  xj  is  the  vector  of  the  generalized  coorditates  of  blade  motion; 

G is  the  matrix  of  the  inertial  coefficients; 

A is  the  matrix  of  rigidity; 

C is  the  matrix  of  the  centrifugal  inertial  coefficients; 

D(t)  is  the  matrix  of  aerodynamic  damping; 

B(t)  is  the  matrix  of  aerodynamic  rigidity. 

If  1^,  ky,  kz  are  infinitely  great,  which  corresponds  to  the  rigid 
fastening  of  the  cyclic  pitch  control,  then  ; ?■  0 and  system  (5)  devel- 
ops into  a flutter  equation  system  for  an  isolated  blade. 

. 

In  those  cases,  when  the  coefficients  of  rigidity  of  the  cables  of 
pitch  and  lateral  control  differ  greatly  from  each  other  or  the  elastic 
attachment  of  the  nonrotating  part  of  the  cyclic  pitch  control  is  such,  ' 

that  the  center  of  rigidity  of  the  attachment  is  located  relatively  far 
from  the  axis  of  the  rotor,  it  is  necessary  to  solve  the  flutter  problem 
of  the  rotor  as  a whole. 
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In  those  cases,  when  the  oscillations  of  a rotor  with  flutter  have 
the  aspect  of  normal  forms  ( for  example,  In  a hovering  mode  with  identi- 
cal rigidities  ^ and  k^) , it  is  possible  to  reduce  a rotor  flutter  problem 
as  a whole  to  an  examination  of  the  oscillations  of  one  isolated  blade  with 
different  values  of  equivalent  rigidity  of  the  control  system,  which  cor- 
respond to  different  normal  forms  of  rotor  oscillations. 


Taking  into  Account  Blade  Elasticity  and  the  Pliability  of  the 
Control  Cables  of  Cyclic  Pitch  Control 

The  method  of  analyzing  the  flutter  in  an  isolated  blade  was  discussed 
in  work  [7],  in  which  employing  the  Bubnov-Galerkin  Method  the  separate 
inherent  forms  of  flexural  and  torsional  oscillations  of  a nonrotating 
blade  as  an  absolutely  rigid  body  were  taken  as  natural  forms  of  flexural 
and  torsional  oscillations.  A calculational  method  was  proposed  in  this  work, 
which  made  it  possible  to  take  blade  elasticity  into  account,  which  was 
attained  by  introducing  additional  flexural  forms  of  oscillations,  obtained 
beforehand  from  frequency  calculation.  In  this  case,  in  order  that  it  would 
be  possible  to  use  the  separate  natural  forms  of  the  oscillations  of  a 
beam  a substitution  of  the  variables  of  the  following  form  is  Introduced 

f (r,  0 — *('.  (6) 

which  made  it  possible  to  separate  the  boundary  conditions  in  the  blade 
root. 

<£**%*-* 


In  this  case,  it  is  assumed  that 
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since  the  moment  of  the  external  forces  relative  to  the  flapping  hinge 
is  much  greater  than  the  moment  of  reaction  in  the  blade  setting. 

In  analyzing  the  flutter  of  an  elastic  blade  employing  the  Bubnov- 
Galerkin  Method  the  zero  and  the  first  natural  form  of  the  flexural  and 
the  first  natural  form  of  the  torsional  oscillations  of  a nonrotating 
beam  in  a vacuum  are  taken.  The  friction  in  the  feathering  hinge  both  for 
an  absolutely  rigid,  as  well  as  for  an  elastic  blade  is  taken  into  consid- 
eration in  the  form  of  additions  to  the  forces  of  aerodynamic  damping. 

As  was  shown  in  a number  of  works  (see,  for  example,  [5,6]),  the 
pliability  of  the  rotor  cyclic  pitch  control,  arising  due  to  the  flexi- 
bility of  the  control  cables,  can  have  a considerable  effect  on  rotor 
flutter,  and  also  the  peculiarities,  as  the  different  rigidity  of  the 
circuits  of  the  pitch  and  lateral  control,  can  even  introduce  qualitative 
changes  in  the  flutter  pattern.  An  attempt  at  analyzing  rotor  flutter  in 
the  flight  mode  taking  into  consideration  the  pliability  of  the  control 
cables  of  the  cyclic  pitch  control  leads  to  great  computational  difficul- 
ties, mainly  because  in  this  case  it  is  necessary  to  investigate  a system 
of  differential  equations  with  periodic  coefficients  of  a higher  order 
than  occurs  for  an  isolated  blade. 

The  lowering  of  the  order  of  this  system  of  equations  by  the  method, 
which  is  employed  in  investigating  ground  resonance,  generally  speaking. 

Is  not  possible.  Work  [5]  substantiated  the  application  of  this  calcula- 
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tional  scheme  by  analyzing  the  accurate  solutions  of  flutter  equations  for 
certain  particular  cases. 

According  to  this  work,  in  equation  system  (5)  the  susbtitution  of 
of  variables  d>f  the  following  form  can  be  carried  out 

*/  +'Xi  COS  Xj  + y,  Sin  (?) 

After  this  substitution  an  equivalent  system  of  differential  equations 
is  obtained  relative  to  the  functions  of  y,  which  contains  periodic  coef- 
ficients of  the  form  sin  KNt  (K  is  a whole  number  and  greater  than  unity). 
Thus,  the  averaging  of  coefficients  with  this  method  leads  to  smaller 
errors,  than  with  the  method  of  averaging  discussed  in  [7].  Thus,  substi- 
tuting (7)  in  (5)  and  orthogonalizlng  the  discrepancy  to  1,  co *Xj,  slnty  , 
in  a set  of  values  of  discrete  variable  j,  we  will  obtain  an  equation  sys- 
tem of  n times  higher  order  than  for  an  isolated  blade  [ n is  the  number 
of  terms  in  representation  (7)]. 

It  is  possible  to  show  [5],  that  as  a result  of  transformation  only 
periodic  coefficients  of  the  type  sinKNt,  cosKNt  can  remain  in  the  new 
equation  system,  where  K is  an  integer,  the  average  of  which  does  not 
affect  the  boundary  of  flutter,  since  the  frequency  of  the  oscillations 
during  flutter  is  close  to  the  frequency  of  the  revolutions.  The  main 
difficulties  in  the  transformation  is  the  process  of  orthogonalizlng  to 
functions  (7).  It  was  proposed  in  work  [5]  to  represent  the  periodic 
coefficients  0(tj)  end  #(xy)  by  trigonometric  polynomials.  With  this  ap- 
proach it  is  possible  to  obtain  a new  equation  system  in  closed  form  only 
for  an  incompressible  gas.  The  consideration  of  compressibility  makes 
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such  a method  of  calculating  analytical  expressions  for  coefficients  of 
matrices  and  B(xJ)  unacceptably  cumbersome. 

However,  in  carrying  out  the  indicated  calculational  method  on  an 
electronic  computer  it  is  proceed  in  the  following  manner.  Let  us  substi- 
tute in  (5) 


Let  us  orthogonallze  the  obtained  discrepancy 


i' 

Y 

; t 
‘ 
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to  1,  cos  •:»,  , sin  -fj  , In  the  set  of  values  of  discrete  variables 

j . We  will  obtain  a system  of  m x n equations  [n  is  the  order  of  the 
matrices  A,  B,  C,  D,  and  m is  the  number  of  terms  in  substitution  (8)] 
with  constant  coefficients. 

Thus,  in  analyzing  an  elastic  blade  for  flutter  taking  cyclic  pitch 
control  into  consideration  it  is  necessary  to  investigate  the  stability 
of  the  solutions  of  a system  of  18  differential  equations  of  the  first 
order  with  constant  coefficients.  If  in  equation  (9)  we  set  Gz  =0, 

Gx  = 0,  Gy  - 0,  then  the  connection  between  the  blades  through  the  cyclic 
pitch  control  disappears  and  system  (9)  now  represents  a system  of  dif- 
ferential equations  of  an  isolated  blade,  but  which  takes  into  considera- 
tion a greater  number  of  harmonics  of  oscillations,  than  in  the  simple 
averaging  of  the  periodic  coefficients  in  equation  (5). 

The  question  of  the  stability  of  the  solutions  of  such  approximate 
equations  reduces  to  an  analysis  of  the  roots  of  a characteristic  equation 
of  this  system.  The  order  of  the  system  of  such  equations,  and  consequent- 
ly, also  the  accuracy  of  the  calculation  of  the  boundary  of  flutter  depends 
on  how  many  forms  of  natural  oscillations  we  take  for  the  isolated  blade 
and  how  many  terms  are  taken  in  substitution  (7). 

An  analogous  procedure  of  approximate  Investigation  of  the  stability 
of  a system  of  differential  equations  with  periodic  coefficients  was  used 
in  work  [9]  in  investigating  the  ground  resonance  of  a helicopter  rotor. 

The  effect  on  the  stability  of  the  motion  and  the  flutter  of  a hell- 
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copter  rotor  of  such  factors,  as  blade  elasticity  and  compressibility  of 
the  flow  incident  on  the  blade  was  investigated  in  this  work,  in  contrast 
to  the  works  indicated  above.  The  maximum  permissible  order  in  such  a 
method  of  the  investigation  of  stability,  achieved  on  an  BESM-3  computer, 
is  12,  i.e.,  nm  12,  where  n is  the  number  of  natural  forms  of  oscilla- 
tions of  an  isolated  blade  (m  = 2 or  3),  and  m is  the  number  of  terms  in 
substitution  (7). 

Results  of  the  Analysis 

Fig.  1-7  show  the  results  of  the  flutter  analysis  of  an  isolated 
blade  and  a rotor  in  accordance  with  a program  compiled  for  the  BESM-3 
computer. 


The  main  parameter,  which  determines  the  margin  of  safety  from  flutter, 
is  the  margin  with  respect  to  parameter  . The  variation  in  the  effec- 
tive centering  of  the  blade  to  the  assigned  was  attained  by  loading  the 
trailing  and  leading  edge  with  a load  corresponding  in  magnitude.  Fig.  1 
shows  the  results  of  the  flutter  analysis  of  a dynamically  similar  model 
of  a helicopter  rotor  blade  of  large  diameter.  The  blade,  having  considerr 
able  elasticity  in  the  plane  of  the  flap  or  stroke  was  tested  for  flutter 
with  different  antiflutter  balancers,  which  were  located  on  the  leading 
edge  at  different  sites  over  the  radius. 

It  was  ascertained  in  the  experiment,  that  the  variation  in  revolutions 
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depends  not  only  on  the  Increase  in  the  effective  centering,  but  also  on 
the  site  at  which  the  antiflutter  balancer  is  located.  The  flutter  analysis 
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of  the  model,  carried  out  taking  three  degrees  of  freedom  into  account, 
confirmed  this  dependence.  The  margin  with  respect  to  aerodynamic  cen- 
tering for  an  elastic  blade  is  equal  to  (where  ; is  the  increase 

I , 

in  effective  centering  for  a semi-rigid  blade  in  the  i-th  cross  section; 
is  a coefficient,  which  depends  on  r^  and  is  approximately  equal  to 

-2 

r^) . For  a semi-rigid  blade  in  accordance  with  the  calculation,  the  criti 
cal  revolutions  depend  only  on  the  the  magnitude  of  the  antiflutter  bal- 
ancer. Fig.  2 gives  the  results  of  the  calculation  of  the  critical  revo- 
lutions for  an  elastic  blade  depending  on  flying  speed.  It  is  evident, 
that  with  an  increase  in  M number  from  zero  to  0.25,  the  critical  revo- 
lutions of  the  blade  fall  off  by  approximately  107..  Fig.  3 shows  the 
change  in  the  decrements  for  elastic  and  rigid  blades  during  transit  of 
the  flutter  zone  boundary.  As  is  evident  from  the  graphs,  the  increase 
in  the  oscillations  in  a blade,  having  considerable  elasticity,  during 
transit  of  the  flutter  zone  boundary,  is  an  order  greater,  than  for  an 
analogous  rigid  blade. 
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Key:  (1)  Calculation;  (2)  Exper- 
iment; (3)  Elastic  blade. 


Fig.  2 

Key:  (1)  Incompressible  gas;  (2)  2 deg 
rees  of  freedom;  (3)  Compressible  gas; 
(4)  3 degrees  of  freedom;  (5)  Compres- 
sible gas;  (6)  2 degrees  of  freedom. 
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Fig.  3 


Fig.  4 shows  the  results  of  the 


flutter  analyses  obtained  by  different 
approximate  methods,  which  take  into 
account  the  different  number  of  harmon- 
ics of  periodic  coefficients  and  gives 
a comparison  of  the  calculations  with 


an  accurate  solution.  It  is  evident 


from  the  results  of  the  calculation, 
that  "simple"  averaging  of  the  periodic 
coefficients*  makes  it  possible  with 
sufficient  accuracy  for  practical  pur- 
poses to  estimate  the  flutter  boundary 
to  a flying  speed  of  Al?^,  0.25. "Three- 
term"  averaging  after  the  introduction 
of  substitution  (7)  makes  it  possible  to 


Key:  (1)  Elastic  blade;  (2)  Com-  . , _ , . 

pressible  gas;  (3)  Rigid  blade;  rather  aCCurately  «rry  out  an  estimation 


(4'  Incompressible  gas. 


of  critical  centering  to  a flight  mode 


of  ft  ■ 0.7—  0.8. 

It  is  also  eveident  from  the  figures  shown,  that  the  critical  center- 
ing, obtained  by  taking  the  periodically  changing  coefficients  into  con- 
consideration,  decreases  more  rapidly  than  the  critical  centering  obtained 
during  calculation  with  averaging  of  the  periodic  coefficients.  For  each 
averaging  method  there  is  a maximum  value  of  the  number  , for  which 
we  will  employ  this  calculational  method.  The  greater  the  number  of  terms 


* 121  is  substitution  of  x ■ y^  (simple  averaging),  122  is  substitution 
of  x • y^  + cos  y*  , 123  is  the  substitution  of  x ■ y + y ^ cos  ^ + 


+ y2  sin  If. 
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in  substitution  (7),  then  for  the  greater  it  is  possible  to  carry  out 
the  calculation.  However,  beginning  with  0.9  it  is  practically 

impossible  to  carry  out  the  calculation  with  the  number  of  terms  in 
substitution  (7)  m > 4,  since  here  computer  error  begins  to  be  ex- 
pressed. The  same  pattern  is  also  observed  in  analyzing  rotot  flutter. 


Mjtjttpttttttt  states* 


Key:  (1)  Isolated  blade;  (2)  Flut 
ter;  (3)  Accurate  solution. 


Kjeteptltttt*  atttcJB* 
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Fig.  5 


Key:  (1)  Isolated  blade;  (2)  Flut- 
ter; (3)  Accurate  solution 


Fig.  6 

Key:  (1)  Rotor;  (2)  Flutter. 


Maam 


Key:  (1)  Rotor;  (2)  Compressible  flow; 
(3)  Incompressible  flow. 


The  practical  range  of  numbers  , in  which  it  is  possible  to  carry 


out  helicopter  rotor  flutter  analysis  is  — 0 — 0.6 
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Fig.  7 presents  the  results  of  the  flutter  analysis  of  a helicopter 
rotor,  which  has  different  control  cable  rigidity.  The  caaculation  was 
carried  out  observing  the  condition  that  the  total  rigidity  of  the  circuits 
of  pitch  and  lateral  control  remained  constant.  Changes  in  the  margin 
with  respect  to  parameter  Ak<.  were  obtained  for  compressible  and  incom- 
pressible gases.  It  is  evident  from  the  figure,  that  calculation  under 
the  assumption  of  incompressibility  qualitatively  reflects  the  pattern 
of  the  variation  in  margins  depending  on  the  relationship  of  the  rigidities 
of  the  control  cables  for  different  flight  speeds.  A sharp  change  in  the 

increase  in  margin  in  the  region  k /k  ■ 1 for  flight  speeds  Af*  + 0 is 

x y 

accompanied  by  a variation  in  flutter  frequency 

Conclusions 

1.  For  a blade,  having  considerable  elasticity  in  the  flapping  plane, 
the  effectiveness  of  an  antiflutter  balancer  depends  not  only  on  its 
magnitude,  but  also  on  the  site  along  the  span,  at  which  it  is  located. 

2.  The  oscillations  of  an  elastic  blade  during  transit  of  the  flutter 
bbundary  increase  by  an  order  more  rapidly,  than  for  an  absolutely  rigid 
blade. 

3.  The  method  of  simple  averaging  of  the  periodic  coefficients  of 


flutter  equations  when  0.25  leads  to  an  exaggeration  of  the  actual 
margins  with  respect  to  the  centering,  which  ensures  safety  from  flutter 
in  flight. 
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